Eigenmodes and thermodynamics of a Coulomb chain in a harmonic potential 



O 

o 

(N 

X) 
D 

oo 



Giovanna Morigi^ and Shmuel Fishman^ 
^ Department of Quantum Physics, University of Ulm, D-89069 Ulm, Germany 
^ Department of Physics, Technion, 32000 Haifa, Israel 
(Dated: February 2, 2008) 

The density of ions trapped in a harmonic potential in one dimension is not uniform. Consequently 
the eigenmodes are not phonons. We calculate the long wavelength modes in the continuum limit, 
and evaluate the density of states in the short wavelength limit for chains of A'^ ^ 1 ions. Remarkably, 
the results that are found analytically in the thermodynamic limit provide a good estimate of the 
spectrum of excitations of small chains down to few tens of ions. The spectra are used to compute the 
thermodynamic functions of the chain. Deviations from extensivity of the thermodynamic quantities 
are found. An analytic expression for the critical transverse frequency determining the stability of 
a linear chain is derived. 
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Cold atomic gases are systems allowing an extremely 
high degree of control, opening new directions and chal- 
lenges both for experiments and theory. The realization 
of low dimensional ultracold gases enables one to design 
and investigate, among others, systems with peculiar ex- 
citations spectra and thermodynamic properties In 
this context. Coulomb crystals constitute a unique form 
of an ordered structure, formed by cold ions in a confin- 
ing potential, that balances the Coulomb repulsion Q. 
The ions vibrate about fixed positions in analogy to the 
situation in an ordinary solid, while the interparticle dis- 
tance is usually of the order of several micrometers, con- 
stituting an extremely rarefied type of condensed mat- 
ter 0, 0, S 0- Variation of the potential permits to 
control the crystal shape as well as the number of ions, 
allowing to explore crystals of very different sizes, thus 
offering the opportunity of studying the transition from 
few particles to mesoscopic systems dynamics Be- 
sides, these structures provide promising applications for 
spectroscopy , frequenc y st andards 9] , study and con- 
trol of chemical reactions [lOj , and quantum information 
processors [lll[l^ll3| . 

In this letter we investigate the dynamics of Coulomb 
chains, i.e. one-dimensional structures obtained by means 
of strong transverse confinement and that usually consist 
of dozens of ions localized along the trap axis [J, |5j . The 
equilibrium charge distribution for such chains is not uni- 
form I14LIl'tI |. This is in contrast to the three-dimensional 
ordering, where the density is uniform and the eigen- 
modes are phonons. In this letter, the chain excitations 
are calculated. These are fundamentally different from 
the phonons in solids because of the non-uniformity of 
the ion distribution and the long range of the Coulomb 
interaction. We derive some thermodynamic functions 
of the chain, and define a specific thermodynamic limit. 
Interestingly, the linear ion chain differs from systems 
that are traditionally studied in the framework of sta- 
tistical physics, since the thermodynamic quantities are 
not extensive. This behaviour is due to the strong cor- 
relation and to the crystal dimensionality. It manifests, 
for instance, in the dependence of the specific heat on 



the number of ions, as we show below. The detailed an- 
alytical derivation of the results presented here will be 
published elsewhere [l3 |. 

We consider TV ions of mass m and charge Q, which 
are confined by a harmonic trap of cylindrical symmetry 
and are crystallized along the trap axis at the positions 
{x"^^ ,0,0), at which the harmonic force and the Coulomb 
force, exerted by the other ions, balance. At sufhciently 
low temperatures the vibrations around these points can 
be considered harmonic and described by 
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where qi = xi — x\ and yi are the displacements in the 
axial and transverse directions (the equations for Zi are 
similar to the ones for y^), while the harmonic confine- 
ment is characterized by the axial and transverse frequen- 
cies V, vt, respectively (yt ^ v)- The coupling matrix is 
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with Wi — qi,yi,Zi. Equations H1I2|I describe a system 
of coupled oscillators, with long range interaction and 
position-dependent coupling strength. In what follows 
we calculate the eigenmodes of this system. For this we 
assume 'Wi{t) = J e"^*?i)i(w)da'/27r. To simplify notations 
we replace Wi{uj) by Wi. This results in equations for the 
eigenmodes of frequency oj that are similar to H1I2|I , but 
with Wi replaced by —uP'Wi. The problem reduces to the 
eigenvalue equation 
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where n = 1,...,N labels the mode. The axial and 
transverse eigenfrequencies are given by the relations 
Lui'^ = v'^ + Xn and lu:^'^ = — A„/2. Since A„ > 0, the 
axial collective excitations are of higher frequencies than 
the trap frequency v, while the transverse excitations are 
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of lower frequencies than vt. In general, Eqs. H1I2|I are in- 
variant under reflection with respect to the center of the 
trap, which coincides with the origin of the axes. Hence, 
the normal modes of the chain are symmetric (even) or 
antisymmetric (odd) under reflection with respect to the 
center, Wi = zLwN-i £3- It can be easily verified that 
the center-of-mass motion is an even eigenmode, w^^\ of 
the Eqs. ^ and (O at frequency v, vt, respectively. Re- 
markably, also the first odd eigenmodes can be exactly 
determined. They are w^"^^ oc x^^"^ at the eigenfrequencies 
y/Siy for the axial, ^Jvl^^iP' for the transverse excitation, 
independent of the number of ions, as can be found by 
substitution of this ansatz in Eq. H3I4|I . combined with 
the equilibrium condition 'iG*!. 

In order to solve Eq. we assume ^ 1. In 
this limit, we may consider the intcrparticle spacing 
ah(xi) = x\^i — xf^ a smooth function of the posi- 
tion, which is inversely proportional to the density of 
ions per unit length nL{xi) — l/aL{xi). The density 
ulIx) is determined assuming a uniform distribution of 
charges inside an oblate ellipsoid of axial length 2L, and 
takes the form nL{x) — SN/AL (l ~ x^ / L^) , which is 
defined for |a;| < L, while the length L is found by 
minimizing the energy of the crystal, and fulfills the 
relation L{N)^ = 3 (Q"^ / mv"^) N log N Q. These ap- 
proximations provide a good estimate of the ions dis- 
tribution in the center of the chain for N sufficiently 
large. In fact, ^^(x) and L{N) used here are the lead- 
ing terms of an expansion in powers of 1/ log N |l4l| . 
For A'^ ^ 1 one can approximate the chain by a con- 
tinuous distribution of charges for the evaluation of the 
long-wavelength modes. In this limit, xf^"^ is a contin- 
uous variable in the interval (—L,L), the displacement 
Wi — 'w{x^) is a continuous function denoted by w(a:), 
and the sum in (PJ is approximated by an integral, where 
the density til (a;) appears as a weight function. In par- 
ticular, the eigenmodes fulfill the orthogonality relation 
J^^dxnL{x)w^"^\x)*w^"^\x) = Sn,m- Usiug an integra- 
tion by parts the sum © can be approximated by 

I[x, w{x)] « -Xo[(l - eW'ix) - ^w'{x)] (5) 

where ^ = x/L and Iq = 3NQ^/{2mL^) « 
The approximation is valid at leading order in log TV 
and sufficiently far away from the edges of the chain. 
The eigenf unctions of Eq. (jSJ are the Jacobi polynomials 
Pe'HO with the eigenvalues A^'^'^ = -£{£ + 3)iy^/2 
with £ = 0,1, . . . and n = £ + 1. Substituting this result 
into the eigenvalue equations (|3I4|I . one finds 

J- = ^^n{n +l)/2. (6) 

Analogously, the eigenfrequencies of the transverse 
modes are = \J v} ~v'^{n- \){n + 2)/4. It 

should be noted that the Jacobi polynomials at f = 0, 1 
are exact solutions of Eqs. H1I2|) . In general, they are 
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FIG. 1: (a) Transverse and (b) axial spectrum. The eigenfre- 
quencies are in units of v. Here, N = 1000 and Vt — AOOu. 
Solid line: numerical solution of Eqs. 11131 . Dashed line: 
(a) Lu^ J"", (b) wil ■'"": These curves have been truncated, as 
they do not correctly reproduce the short wavelength eigen- 
modes. Dashed-dotted line: spectra obtained from the den- 
sity of states 0. 

the correct eigenmodes of Eq. I^J in the asymptotic limit 
N oo. For finite N a perturbative expansion in the 
parameter 1/ log TV is required, resulting in a very slow 
convergence. Nevertheless, the behaviour A„ (x i'^ is ex- 
act in the continuum limit 0,^3- Figure ^ presents the 
comparison between the spectrum of eigenfrequencies, 
calculated numerically for 1000 ions, and the solutions 
Lofi Figure 12 shows the relative deviation of 

the frequencies ^ from the numerical results for chains 
of different number of ions. From Fig. |2Ia) one sees that 
Eq. ® gives already a good estimate of the eigenfrequen- 
cies for a chain of 10 ions. Its prediction is valid for the 
axial low- and the transverse high-frequencies, and it im- 
proves slowly as N increases, as expected from the slow 
convergence of the 1 / log N expansion. 

The short wavelength eigenmodes are characterized by 
relatively large displacements of the ions at the chain 
center, while the ions at the edges nearly do not move. 
In fact, the intcrparticle distance is minimal around the 
center, and for TV ^ 1 it is of order 1/TV, while it is 
significantly larger at the edges. Hence, a wave cannot 
propagate to the edges, where the intcrparticle spacing 
is larger than the wavelength. Moreover, in the center 
and for the short-wavelength excitations we expect that 
the relevant contributions to the force on an ion originate 
from the neighbouring charges, as nearby groups of ions 
move in opposite directions, resulting in forces that mu- 
tually cancel. By this hypothesis we keep in Q only the 
nearest-neighbour interaction and apply the formalism 
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FIG. 2: (a) SLOn/u)n and (b) 5a;i/Aa;i as a function of n, 
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cies Ljjl are obtained by solving numerically From top 

to bottom: = 10, 50, 200, 500, 1000. 
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FIG. 3: acr(A'') as a function of the number of ions. The solid 
curve gives the result according to . The dashed line is the 
curve 2.53 N~^'''^ , that fits the values obtained by molecular 
dynamics simulations |23.|2^. 



developed by Dyson |20| to derive the density of states 
as a function of the squared frequencies A. This is given 
by 

/(A) 



D{X)=8/{nNX) dxTiLix) + 4:A{x)/ X (7) 



where A{x) = 2Q'^nL{x)^ /m is the nearest-neighbour 
coupling constant for slowly-varying inter par ticle dis- 
tance, and /(A) = - (A/4A(0))i/3 Here, 
A(0) « 9iy'^N'^/32logN at leading order in logiV. The 
spectrum obtained from Eq. lO is shown in Fig. ^ The 
deviations from the numerical result are due to the as- 
sumption of nearest-neighbour coupling, and are small 
in the short-wavelength part of the spectrum, show- 
ing that Eq. (Q provides a good approximation in this 
regime. In particular, it provides a good estimate of 
the eigenvalue Ajv, which determines the maximal axial- 
and minimal transverse-frequencies through the relations 
wJlax = V*^^ + Aat and tu^i^ = ^v} - Ajv/2. This eigen- 
value is found from /(A) = 0, as for A > Aat the function 
_D(A) vanishes, and Aat = 9z^^iV^/(81og A^) at leading or- 
der in logA^. The minimal transverse frequency is 
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It can vanish for certain values of vt and N . We identify 
the critical value ft,cr = \/ Atv/S, such that for vt < vt,ci 
the linear chain is unstable. The critical aspect ratio 
between the trap frequencies Ofcr — ^^1^1^^ takes the 
value 
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which fixes the condition on vt for the chain stability ac- 
cording to the inequahty > a^^v^ . The result Q is in 
agreement with the analytical estimate in |2l| , which was 
obtained under different requirements. Figure 13 shows 
that it also approximately reproduces the curve obtained 
by fitting points calculated by molecular dynamics sim- 
ulation |22| . which has been experimentally verified for 



chains of tens of ions [23. It was suggested that this 
instability, resulting in a structural transition from a lin- 
ear to a zig-zag equilibrium confi gur ation 0, l22l | , can be 
treated as a phase transition |2^ 12,'^ . 

Using these results we now discuss the thermodynam- 
ics of a linear crystal of iV ions, whose dynamics are 
described by SA'^ harmonic oscillators of frequencies wl, 
cj^, in the regime where the axial frequencies are much 
smaller than the transverse excitations, i.e. wjiiax ^ ^min- 
The eigenmodes are quantized using standard proce- 
dures 'l7'|. The chain is assumed to be in thermal equi- 
librium with a bath at temperature T, which is suf- 
ficiently low so that only the axial modes are excited 
(ksT <C fiuj^^^. Hence, the dynamics of the system is 
one-dimensional. The thermodynamic limit for this kind 
of system can be defined by assuming constant interpar- 
ticle spacing (thus constant linear density) at the cen- 
ter of the chain as A" increases, namely requiring ai(0) 
to be constant, in analogy with the definition for cold 
gases in traps Since 0^(0) oc (-^log N /vN)'^^^ , this 
requirement corresponds to a vanishing axial frequency 
as A^ 00, according to v ^ \/\ogN /N. Using this 
definition, wLac and w^jn take finite values in the ther- 
modynamic limit, since they depend on A^ and v only 
through the combination fN/y/logN. We identify the 
thermodynamic variables with T, and i/, whose vari- 
ation results in a variation in the length of the chain j3| , 
and take constant First we study the heat capac- 
ity Ca = dUth/dT\^,N, where the thermal energy Uth is 
the average energy of the excited state relative to the 
ground state energy. At high temperatures Ca exhibits 
the Dulong-Petit law Ca = Nks, while at low tempera- 
tures it takes the form 
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where xq = [3h,v and we have used the density of states 
derived from ©. For ksT :§> Tiv we can set ~ and 
obtain Cq oc T . Hence, in this temperature regime the 
heat capacity is linear in T, like for a one-dimensional 
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FIG. 4: Coefficient of tliermal expansion qt, in arbitrary 
units, as a function of T, in units of Tiuj\i/kB- Here, A'^ = 1000 
and ojmax '^min- This would correspond to an experiment 
with V — 2-K X 1 kHz and Ut = 2-k x ?> MHz. 

Debye crystal j23]. Remarkably, for these temperatures 
Co is inversely proportional to v. Hence, the specific 
heat per particle Ca — Ca/N behaves like Ca ~ \/Nv and 
vanishes as Ca ~ 1 / Vlog N in the thermodynamic limit 
defined above. Thus, at low temperatures it depends on 
the number of ions, thereby manifesting a deviation from 
extensivity of the system's behaviour. 

The coefficient of thermal expansion ar also exhibits a 
remarkable behaviour as a function of the temperature. 
It is related to Ca and to the isothermal compressibility 
Kt according to the relation ax = ZHirCa/'^L [i^. In 
the regime of thermal stability, the isothermal compress- 
ibility does not vanish and is practically independent of 
the temperature, the pressure being dominated by the 
zero-temperature component. In particular, it scales as 
Kt ~ 1/logiV. Hence, the temperature dependence of 
ax is solely determined by the heat capacity Ca- At low 
temperatures Ca l/i', and as the thermodynamic limit 
is approached ar ''^ (log A^)~'^/^. At higher tempera- 
tures, when the heat capacity manifests the Dulong-Petit 
behaviour Ca = Nks, the coefficient of thermal expan- 
sion vanishes like ar ~ 1/ log N. For large but finite 
number of ions ax does not vanish, in contrast to the 
behaviour found in a uniform harmonic solid j^^. The 
coefficient of thermal expansion ut is plotted in Fig. 4 
as a function of T for a chain of 1000 ions. 

Finally, we remark on the regime of thermal stability, 
which holds when the thermal energy is considerably 
smaller than the equilibrium energy, or, equivalently, the 
displacements are much smaller than the respective inter- 
particle distances. The condition for thermal stability is 
\ogN/aL{0) ^ ksT. Outside this regime, thermal ex- 
citations may cause structural transitions, accompanied 
by a critical behaviour of the thermodynamic functions. 

In conclusion, we studied the low temperature dynam- 
ics of a strongly-correlated and one-dimensional crystal 
inside a harmonic potential. We investigated how the 
dynamics of the system is affected by the dimensionality 
and the long-range correlations and how this is manifest 
in the thermodynamics functions. In future works we will 
explore the behaviour of these functions at the instabil- 



ity points, in order to characterize the crystal structural 
transitions and their relations to the standard theory of 
phase transitions |2^ . 

The excitations spectra we presented are of interest for 
spectroscopy and quantum information with ion traps, 
where quantum logic is imp lemented with the lowest 
collective modes j^, 0|. This work may be rele- 
vant to studies of Coulomb systems like ion cry stals in 
storage rings and cold neutral plasmas which 
at sufficiently low temperatures are expected to crystal- 
lize , and contributes to the on-going research on low- 
dimensional cold gases 0]. 
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